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jy^ . Coalescing black-hole binaries are expected to be the strongest sources of gravitational waves for 

ground-based interferometers as well as the space-based interferometer LISA. Recent progress 
in numerical relativity now makes it possible to calculate the waveforms from the strong-field 
dynamical merger and is revolutionizing our understanding of these systems. We review these 
dramatic developments, emphasizing applications to issues in gravitational wave observations. 
These new capabilities also make possible accurate calculations of the recoil or kick imparted 
i' to the final remnant black hole when the merging components have unequal masses, or unequal 

JL ' or unaligned spins. We highlight recent work in this area, focusing on results of interest to 

astrophysics. 
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1. Introduction 

Gravitational wave astronomy will open a new observational window on the universe. 
' Since large masses concentrated in small volumes and moving at high velocities generate 

the strongest and therefore most readily detectable waves, the final coalescence of black- 
hole binaries is expected to be one of the strongest sources. The opening of the full 
electromagnetic spectrum to astronomical observation during the last century greatly 
expanded our understanding of the cosmos. In this new century, observations across the 
gravitational wave spectrum will provide a wealth of new knowledge, including accurate 
measurements of binary black-hole masses and spins. 

The high frequency part of the gravitational wave spectrum, ~ 10Hz < / < 10 3 Hz, 
is being opened today through the pioneering efforts of first-generation ground-based 
interferometers such as LIGO, currently operating at design sensitivity. Such instru- 
ments can detect gravitational waves from coalescing stellar-mass (M < 1O 2 M ) and 
intermediate-mass (10 2 M© < M < 10 3 A1©) black-hole binaries. While detections from 
this first generation of detectors are likely to be rare, the advanced LIGO (adLIGO) 
upgrade may detect the coalescence of several stellar-mass and tens of intermediate-mass 
black-hole binaries per year. Other high-frequency sources include binary neutron star 
coalescences, supernovae, and rotating neutron stars. 

The low frequency gravitational wave window, 3 x lCU 5 Hz < / < 1Hz, is especially rich 
in astrophysical sources and will be opened by the space-based LISA detector, currently 
in the formulation stage. LISA will be sensitive to the coalescence of massive black-hole 
binaries with total masses in the range 10 4 M Q < M < 10 7 M Q to large redshifts z > 10 at 
relatively high signal-to-noise ratios (SNRs), and may detect 10 or more such events per 
year. Using such observations, the black-hole masses, spins and luminosi ty distances can 
be de termined to very good precision, with errors < 1% in some cases (jLang fc Hughes! 
2006). In addition, LISA will detect gravitational waves from the inspiral of compact stars 



into central massive black holes out to z ~ 1, as well as tens of thousands of compact 
binaries in the Galaxy. 
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The actual merger of two comparable-mass black holes that plunge together and form 
a common event horizon takes place in the strong-field dynamical regime of general rel- 
ativity. For many years, we were unable to calculate the expected waveforms from these 
very energetic events due to severe problems with the large scale computer codes needed 
to simulate the mergers. Recently, however, a series of stunning breakthroughs has oc- 
curred in numerical relativity, resulting in stable, robust, and accurate simulations of 
black-hole mergers as well as applications to astrophysics. In Sec. [5] we review these de- 
velopments and present examples of the resulting gravitational waveforms. Applications 
of these signals to issues in gravitational wave observations are discussed in Sec. [3] When 
the merging black holes have unequal masses, or unequal or unaligned spins, the final 
remnant black hole suffers a recoil; recent progress in calculating these "kicks" and their 
applications to astrophysics are presented in Sec. [U We conclude with a summary in 
Sec. El 



2. Calculating black-hole binary coalescence 

The final coalescence of a black-hole binary is driven by gravitational wave emission, 
and pr oceeds in three stages: an adiabatic inspiral, a dynamical merger, and a final ring- 
down ([Flanagan k. Hughedll998h . During the inspiral, the black holes are well-separated 
and can be approximated as point particles. The black holes spiral together on quasi- 
circular trajectories, and the resulting gravitational waveforms are chirps, i.e. sinusoids 
that increase in both frequency and amplitude as the black holes get closer together. The 
inspiral can be treated analytically using the post-Newtonian (PN) ap proach, which is 
an expansion in v/c, where v is the characteristic orbital velocity (see iBlanchet ( 2006h 
for a review of PN results). The inspiral is followed by a dynamical merger in which 
the black holes plunge together to form a highly distorted single black hole, producing 
a powerful burst of gravitational radiation. Since the merger stage occurs in the regime 
of very strong gravity, a full understanding of this process requires numerical relativity 
simulations of the Einstein equations. After merger, the remnant black hole then settles 
down, evolving towards a quiescent Kerr state by shedding its nonaxisymmetric modes as 
gravitational waves. The late part of this ringdown stage can be treated analytically using 
black-hole perturbation theory, and the resulting gravitational waveforms are superposi- 
tions of exponentially damped sinusoids of constant frequency ( Leaver 19861 : Echeverrfal 
19891) . 



In numerical relativity, the full set of Einstein's equations are solved on a computer in 
the dynamical, nonlinear regime. This is typically accomplished by slicing 4-D spacetime 



into a stack of 3-D spacelike hypersufaces, each labeled by time t ([Arnowitt et q/.lll962 



Misner et aLlll973 ). The Einstein equations split into two sets. The constraints give a set 



of relationships that must hold on each slice and in particular constrain the initial data 
for a black-hole binary simulation. This data is then propagated forward in time using 
the evolution equations. Four freely-specifiable coordinate, or gauge, conditions give the 
development of the spatial and temporal coordinates during the evolution. 

Simulating the merger of a black-hole binary using numerical relativity has proved to 
be very challenging. The first attempt to evol ve a head-on collision in 2-D axisymmetry 
dates back to 1964 ( Hahn k, Lindquist 19641 ). In the mid-1970s, the head-on collision 
of two equal mass, nonspinning black holes was first simulated succe ssfully, along with 



1976; 



the extraction of s ome information about the gravitational radiation (|Smarr et 

19791) . In the 1990s, fully 3-D numer ical relativity codes were developed 



Smarr 1977, 



and used to evolve gra zing collisions of black holes ( Brugmarml [l999; Br andt et al. 2000l : 



Alcubierre et a/.ll2001l ). However, the codes were plagued by a host of instabilities that 
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caused them to crash before any significant portion of a black-hole binary orbit could be 
evolved. For many years, progress was slow and incremental. 

Recently, a series of dramatic developments has led to major progress in black-hole 
binary simulations across a broad front. The firs t complete orbit of a black-hole bi- 
nary was achieved in 2004 (jBrugmann et all 120041 ). This was followed by the first full 



simul ation of a black -hole binary through an orbit, plunge, merger and ringdown in 
2005 (|Pretoriusll2005h . In late 2005, the development of new coordinate conditions pro- 
duced a breakthrough in t he ability to carry out accurate and stable long-term evolu- 



tions of black-hole binaries (jCampanelli et aZ.l l2006a: B aker et a l. 2006c; van Met er et al 



l2006h . These novel but simple "moving puncture" techniques proved highly effective 
They were quickly adopted by a broad segment of the numerical relativity community, 
leading to stunning advances in black-hole binary modeling, starting with evolutions 
of equal mass, nonsp inning black holes and m o ving quickly t o inclu d e unequal masses 
and spi n s; (see, e.g. , Campanelli et al. ( 20066 ): IBaker et al\ (20066 ): Campanelli et al 
(12006 J: [Gonzalez et all (1200761): iBaker et all (1200761): iHerrmann et all (120061 l2007fh 
ICampanelli et all (I2007af6lh iKoppitz et all (|2007l ): iGonzalez et all (|20o'7al ): IBaker et al 
(|2007a[ ): lTichv fc Marronettil (|2007i) ~ 

The most rapid advances in modeling black-hole binary coalescences cover the previ- 
ously least understood part of the gravitational waveform, i.e. the final few cycles of ra- 
diation generated from near the "innermost stable circular orbit" (ISCO) and afterward, 
which we call the "merger-ringdown" . There is already considerable progress toward a 
full understanding of this important "burst" portion of the waveform, through which the 
frequency sweeps by a factor of ~ 3 up to ringdown and during which the gravitational 
wave luminosity is ~ 10 23 L©, more than the luminosity of the combined starlight in the 
visible universe. 

A particularly significant development was the demonstration of init ial- d at a- indep en- 



dence of merger-ringdown waveforms for equal-mass, nonspinning black holes ( B aker et 



200661 ). as summarized in Fig. [TJ Results from four runs with successively larger initial 



separations are shown; the waveforms have been aligned so that the moment of peak 
radiation amplitude in each simulation occurs at time t — 0. Here we show the gravita- 
tional wave strain from the dominant (I = 2, m = 2) mode; this represents an observation 
made on the equatorial plane of the system, where only a single polarization component 
contributes to the measured strain. The upper panel of Fig. [1] shows the full simulation 
waveforms, while the lower panel focuses on the final merger-ringdown burst. Note that 
here and elsewhere in this paper, we use geometrical units, G = c = 1, to measure time, 
distance and mass in the same units. In particular, one solar mass Mq is equivalent 
to ~ 5 x 10 sec, or ~ 1.5kmQ In Fig. [lj our timescale is the final mass m/, of the 
post-merger hole; this will be less than the initial total binary mass M because of grav- 
itational radiation. In the shortest run (solid line), the black holes are placed on initial 
orbits close to the ISCO and undergo a brief plunge followed by a merger and ring- 
down (jBaker et aZ.ll2006dl . At the next-larges t initial separation (dashed line), the black 
holes complete ~ 1.8 orbits before merging ( Baker et aLll20066h . The waveforms from 
the two runs with successively larger initial separations (dot-dashed and dotted lines, 
respectively) then lock on to the merger-ringdown part of this shorter (dashed) run. In 
fact, the waveforms from the three longer simulations show very strong agreement for 
t > —50 m/, with differences among these waveforms < 1% in this regime. 

Today, there is consensus that the merger of equal-mass, nonspinning black holes pro- 



f Since the simulation results scale with the masses of the black holes, they are equally 
applicable to LISA and ground-based detectors. 
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Figure 1 Gravitational waveforms are shown for four simulations of equal mass nonspin- 
ning black-hole binaries with successively wider initial separations. Here, di is the initial 
coordinate separation of the black- hole centers ( "punctures" ) . All waveforms have been 
aligned so that the moment of maximum radiation amplitude occurs at t = (phase is 
also set to zero at t = 0). The upper panel shows the full extent of all the runs, while 
the lower panel shows the final merger-ringdown burst. 



duces a final Kerr black hole with spin a/m/ ~ 0.7, and that the amount of energy 
radiated in the form of gravitational waves, starting with the final few orbits and pro- 
ceeding through the plunge, merger and ringdown is J? ra .ri = M-mf ~ 0.04M (IPretoriui 



20051: Campanelli et all 2006a ; Baker et al. 2006c : Campanelli et al. 2006 \k Baker et al 
20066); see Centrella ~(|2006l) for a review. There is also agreement on the overall sim- 



ple shape of the waveforms shown in Fig. [T] and detailed compariso n of results among 
numerical relativity groups has already begun (jBaker et all 12006 of ) . Parameter-space 
exploration of the merger-ringdown burst, using simulations w ith variou s mass ratios 
(Herrmann et "fllll2006l Baker et al 1l2006o!l : lGonzalez et olll2007fch and spins (jCampanelli et 
200633). 

is now underway. The next step is to push this frontier to increasingly complex 
mass-ratio and spin-orientation combinations, and to establish initial-data independence 
across these parameters. 

Simulations starting in the late-inspiral regime and covering more than a few orbits 



prior t o merger are more challenging than shorter merger-ringdown evolutions (jBaker et al 



200761 ). Such long simulations have stronger requirements for numerical stability and 
place greater demands on computational resources. In addition, better accuracy is needed 
to control the accumulation of phase error; this in turn constrains the numerical error 
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that can be tolerated in the rate of energy loss through gravitational radiation reaction, 
which governs the inspiral. 

So far, the longest simulations have been carried out for equal-mass, nonspinning black- 
hole binaries, s tarting at relatively wide separation s. The two longest r uns shown in Fig.fT] 
undergo ~ 4.2 ( Baker et aLll20066 ) and ~ 7 orbits ( Baker et q/.ll2007fel). respecti v ely, pri or 
to mergeiQ and demonstrate progress in long simulations (see also Husa et al. ( 2007oh ). 
We expect that astrophysical black-hole binaries in the late inspiral regime will follow 
nearly circular orbits, since any initial eccentricity would have been radiated away early in 
their evolution. However, in the first long black-hole binary simulations, the eccentricity 
of the initial orbits was not very well controlled; in particular, the early part of the 
waveform shown by the dot-dashed line shows evidence of eccentricity which resulted from 
the initial data specification. The dotted cu rve shows more rece n t results, which sta rt 
with very small eccentricity e < 0.01 (see also lPfeiffer et all (|2007t) ; lHusa et all (|2007^ 1. 
Such long waveforms make it possible to compare the results of numerical relativity 
simulat ions of black-hole binar i es with post-Newtonian calculations in the late-inspiral 
regime (jBuonanno et all\2007at iBaker et aLll2006et lHannam et aLll2007j) . demonstrating 
remarkable agreement. 



3. Observing Black-Hole Binary mergers 

The final moments of a black-hole binary merger produce the most intense radiation 
generated through the strongest gravitational dynamics. The resulting waveforms are 
obviously an important part of the observable gravitational wave signature for black- 
hole binary events. Before these recent advances in numerical relativity, however, little 
was confidently known about these signals, so that observational anal yses of gravitationa l 
wave data could only rely on generic "unmodeled burst" techniques (| Abbott et q/.ll2007 ) 
rather than more effective matched-filtering techniques, which require detailed knowledge 
of the expected signals. Consequently the quality of scientific information which could be 
expected from merger observations has, to date, necessarily been significantly discounted 
in planning observational work. 

Just how much more we can learn from observing the final stages of black- hole binary 
mergers with full knowledge of the waveforms, and how best to apply this knowledge in 
gravitational wave data analysis, are still largely unanswered questions, which the gravita- 
tional wave community is only beginning to address. We can get some rough impressions 
of the observational significance of black- hole binary mergers from this new perspective by 
considering the S NR of the full waveforms b ased on matched-filtering analysis of gravita- 
tional wave data (jFlanagan fc Hughesl 1998 ). Estimating the full waveforms by stitching 
together the most accurate PN model with a numerical simulation of the last several 
cycles , we have calculated the SNR for equal-mass nonspinning binaries ( Baker et all 
1200761 ). Fig. [5] shows the contours of SNR as functions of redshift z and total binary 
mass M for the ground-based adLIGO detector, which will be sensitive in the frequency 
range lAHz < / < 10 3 Hz. AdLIGO is a planned upgrade of the initial LIGO detectors 
that will increase the sensitivity by roughly an order of magnitude across the frequency 
band. In addition, adLIGO can be tuned to optimize its sensitivity for different sources. 
To produce Fig. [2] we used the wide-band tuning typically associated with burst sources, 
due to its greater sensitivity a t higher frequenci es, where the merger portion from many 
sources is predicted to occur (|Shoemakerll2006h . This gives an improved SNR for most 



f For these quadrupolar waves, the gravitational wave and orbital frequencies are related by 

/ = 2/ or b . 
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Figure 2 Contours of SNR for observations of black-hole binaries with adLIGO, showing 
mass, redshift, and luminosity-distance dependence. 



black-hole masses compared to tunings that were optimized for only the early inspiral 
portion of the coalescence. In Fig. [3] we plot contours of SNR for LISA observations, 
showing that LISA can observe massive black-hole binaries at high SNR throughout 
the observable universe. LISA is most sensitive to systems with masses in the range 
10 5 M Q ^ M ^ 10 7 M©, coinciding with the mass range in which models of massive 
black -hole binaries predict coalescence within a Hubble ti me dMilosavljevic fc Merritt 
2003) and event rates for LISA of at least several per year (jSesana et aZ.ll2005l ). 

In Fig. [4] we plot a simulated LISA data stream containing a black- hole binary signal 
as well as noise from the interferometer and the white dwarf confusion noise caused 
from the unresolvable number of background white dwarf binaries. This is an example 
of a typical, non-optimal case in order to demonstrate the relative strength that these 
signals will generally have. Here we show two 4 x 10 5 Af© black holes at redshift z = 10, 
inclination i — ir/2, polarization angle tp = 0, ecliptic latitude = 0, and ecliptic 
longitude A = ir. So far we have restricted our consideration to the example of equal- 
mass nonspinning black-hole binary mergers. More detailed considerations of black-hole 
binary observations must take into account the full parameter space of these systems, 
which, from the point of view of general relativity, are described by seven parameters, 
corresponding to the black- hole mass-ratio q = m\/m,2 ^ 1 (recall that the total mass M 
only scales the waveforms), and the spin vector of each hole (plus two more parameters 
if eccentric systems are considered). Much progress has already been made in conducting 
n umerical simulations to f lesh out the basic fe a tures of this parameter space (see, e.g. 
, ICampanelli et all l|2006cO ): iHerrmann et all 42006)). First impressions indicate that 



the waveforms in this parameter space are qualitatively simple, with smooth parameter 
dependencies. This suggests that our knowledge of the mergers can ultimately be applied 
in observations with the aid of smooth analytic waveform models, joining PN treatments 
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of the inspiral signals with analytic approximations to the mergers. Already such a model 
has been proposed, whic h provides high overlaps with the numerical simulations for mass- 
ratios as low as q = 1/4 (jBuonanno et a/.ll20076l). 

Applying the mass scalings from Flanagan fc Hughesl ( 19981 ). we can estimate the effect 
of the mass ratio q on the computed SNRs; specifically, SNR ~ -q 1 / 2 for the inspiral, and 
SNR ~ 77 for the merger and ringdown, where 77 = mi?Ti2/M 2 = q/(l + q) 2 is the 
symmetric mass rati o. For stellar-mass bla ck-hole binaries, the mass ratios are rather 
broadly distributed (jBelczvnski et all 120071) : the rates for such mergers may be low, ~ 
2yr~ 1 for adLIGO, depending on the evolution of the original binary through the common 
envelope phase. For intermediate-mass black- hole binaries, mass ratios in the range 0.1 < 
o < 1 are expected to be the most relevant, with potential rates of ~ 10 per year 
( Fregeau et a/.ll2006 ); we note, however, that these rates are far more uncertain than those 
for stellar-mass black-hole binaries. Investigations of gravitatio nal-wave data an alysis 
based on unequal-mass numerical simulations are now underway: IPan et all (|2007l ) have 
investigated how numerical relativity waveforms may be used to select the best post- 
Newtonian template banks for searches in the LIGO, adLIGO and VIRGO data stream. 

Astrophysical black-hole binaries are also expected to be spinning, which can poten- 
tially raise the SNR, for example if the re is a spin-induced hangup that generate s more 
gravitational wave cycles in the merger ( Campanelli et al. 2006(il ). IVaishnav et all (|2007l) 
have investigated how the addition of spins to equal-mass waveforms affects search tem- 
plates, necessitating the use of higher-order modes for high matches in the LIGO band. 

We also note that, with more accurate computational modeling, simulations of the 
merger can be used together with gravitational wave observations to probe gravity in the 
regime of strong fields. In particular, if the binary masses and spins can be obtained with 




Figure 4 Simulated LISA data stream showing LISA's response to a system of two equal- 
mass black holes. Specifically, the pseudo-A observable is shown for the case of two 
4 x 10 5 M© black holes at redshift z — 10, inclination i = tt/2, polarization angle ip = 0, 
ecliptic latitude (3 = 0, and ecliptic longitude A = n. The pseudo-A, E, and T observables 
are a set of linearly inde pendent combinations of the Michelson unequal arm X, Y, and Z 
observables at TDI 1.5 (jShaddock et a/.ll2003 ). and are therefore useful in data analysis. 
The inset shows a time span of 4 x 10 4 seconds encompassing the time span of the main 
figure, the main figure being an enlargement of the dash-boxed region of the inset. 



good accuracy from observations of the inspiral with LISA (|Lang fc Hughes 2006), the 
merger waveform can be calculated using numerical relativity. With this, a comparison 
can be made between the predictions of general relativity - or any other theory of gravity 
used in a numerical simulation - with observations in the regime of very strong gravity. 



4. Kicks from black-hole binary mergers 

Black-hole binaries that possess an asymmetry, such as having unequal masses or 
possessing unequal or unaligned spins, will emit gravitational radiation asymmetrically. 
This asymmetric gravitational radiation will impart a linear momentum recoil on the 
binary's center of mass. Previous atte mpts to calculate this recoil , or "kick" , an alytically 
in the nonspinning case ( Fitchettlll983h and for cases with spin (|Kidderlll995h included 
very large uncertainties, because the dominant part of the kick occurred in the last 
orbit through the merger, which was very poorly modeled. With the advent of stable, 
accurate numerical relativity codes, the problem of accurately calculating the kick became 
possible. Several groups have tackled this problem . The pr oblem of kicks from an unequal 
mass binary was first investigated in lBaker et al. ( 2006rf|). where recoils betwee n 86 and 
115 km/s for a q — 2/3 mass ratio were simulated. In Ibonzalez et ~al\ ( 20076 ). several 
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simulations were performed and a peak recoil velocity for nonspinning binaries of 175 
km/s was found for a mass ratio of q w 1/2.78. 

A flu rry of activity began w hen several groups began simulating kicks from spin asym- 
metry. [HerrmamTljFaZ ( 2007 ) demonstrated that kicks from spin asymmetry can exceed 
unequal-mass kicks, and that for an equal-mass configuration with one black hole spin- 
ning prograde and the other spinning retrograde with equal magnitude a (a dimensionless 
spin parameter, varying from for a Schwarzschild hole to 1 for a maximally spinning 
Kerr hole), the spin kick scales as 475 a km/s. However, it was soon found that the 
largest kicks are obtained when the spins are antiparallel with each ot her but perpendic 



ular to t he orbital angu l ar mom entum, a configuration suggested by ICampanelli et al 



( 2007a ). Gonzalez et al. ( 2007 a J) were the first to perfo rm this simulation, obtaining a 
kick of 2650 km/s for a — 0.8. ICampanelli et all (|20076l ) showed that the resulting kick 
depends sensitively on the angle between the black-hole spin vectors and the i r veloc - 
ities. They found a simple cosine dependence, indicating that IGonzalez et ali (|2007al ) 
somewhat fortuitously ha ppened upon what is nearly the maximum angle. Furthermore, 
ICampanelli et al. ( 20076 ) projected that for maximal spin, kicks as large as 4000 km/s 
are possible. 

Kicks of 2000 km/s are sufficiently large to kick remnant black holes out of even 
the largest giant elliptical galaxies. Therefore, simulation results this large had to be 
reconciled with the observational reality that we see massive black holes at the centers 
of all the galaxies that we hav e observed above a certain spheroidal mass threshold. 
Schnittman fc Buonannol ( 2007 ) used the EOB formalism to perform a Monte Carlo 



simulation in order to try and predict how typical these large kicks are. For mass ratios 
between q = 1 and q = 1/10 and spins of a, = 0.9 on both holes, they allowed random 
orientations among the spin vectors and the orbital angular momentum. They found 
that 12% of the remnants received a kick exceeding 500 km/s, and 2.7% received a kick 

e xceeding 1000 km/s. 

Schnittman fc Buonannol ( 2007 ) assume no external influence which might tend to 



align the spins with each other and with the orbital angular momentum. However, 
Bogdanovic. Re ynolds & Miller (|2007l ) investigated the role of the Bardeen-Petterson ef- 



fect ( Bardeen fc Pettersonlll975 ) in mitigating the gravitational recoil from spin asymme- 



try. In the absence of gas accretion, the assumption of random orientation in lSchnittman fc B 



(2007) applies. However, if a black-hole binary accretes (1-10)% of its initial mass from 
an accretion disk, the spin of each hole will align with the angular momentum of the disk, 
which in turn is aligned with the orbital angular momentum of the binary. Therefore, 
if the majority of mergers occur in gas-rich environments, then the giant kick results, 
although theoretically interesting, are not astrophysically relevant. There is good reason 
to expect supermassive black holes to accrete substantial mass, particularly at larger red- 
shift. In this scenario, the configuration of greatest interest is spin orientations aligned 
with the orbital angular momentum. 

Fig.[5]and Table [1] present the results of our inv estigation of this class of configurations, 
which were first presented in Baker et al. ( 2007al ). "NE" refers to unequal mass, with the 
corresponding simulations having a mass ratio q = 2/3. "+" means prograde with respect 
to the orbital angular momentum, and "-" means retrograde ("0" means no spin). Fig. [5] 
shows the accumulated kick as a function of time for all of our runs. Of particular note 
is the absence of an "unkick" for the equal-mass case, and the variable size of the unkick 

in the other cases. Also, we observe the expected accelerated merger in the NE and 

NEo- cases due to spin-orbit attraction, and correspondingly the delayed merger in the 
NE + + and NEq_i_ cases due to s pin-orbit repuls i on. U sing our data along with the data 



from lHerrmann et al. ( 2007 ) and Koppitz et ali (|2007l ) (which are also included in Table 
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Table 1. Predicted ver sus computed kick speed. Runs labeled " S0.##" are taken fro m 
Herrmann et wh ile runs labeled "r#" are taken from iKoppitz et a/1 (|2007l ). 

[from iBaker et all (|2007af ) - reproduced with permission of the AAS] 
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[l}, we are able to construct an empirical kick formula, given by: 

v = V {32 q 2 /{l + q) 5 } V(l ~ q) 2 + 2 (1 - q) K cos 6 + K 2 , 



(4.1) 



where K = k {qd\ — 62). The parameter Vo gives the overall scaling of the kick (note that 
the factor in brackets becomes unity for q — 1), while k gives the relative scaling of the 
kick contributions from spin and mass asymmetries. 

Performing a least-squares fit to the data yields Vo — 276 km/s, q — 0.58, and k — 0.85. 
Our formula yields a maximum error of 10.8% for the cases investigated. If, in fact, 
the majority of supermassive black-hole mergers are gas-rich, then this formula may 
predict a significant component o f the kick for the majority of cases of astrophysical 
interest. Campanelli et al. (1200761) subsequently suggested an extension of the formula 



presented in IBaker et al\ (|2007al ) for kicks out of the orbital plane, which would be 
more applicable in gas-poor cases where the Bardeen-Petterson effect doesn't occur, 
such as gas-poor galaxies and the merger of stellar-mass black holes. The achievement of 
simulating the merger is therefore leading us to a fairly comprehensive understanding of 
the recoil velocity parameter space in a strikingly brief period of time. 



5. Summary 

Mergers of comparable-mass black-hole binaries produce intense bursts of gravita- 
tional radiation and are very strong sources for both ground-based interferometers and 
the space-based LISA. Recent progress in numerical relativity now makes it possible to 




t/M 



Figure 5 Aggregated kicks from all of our runs. The merger time for each binary matches 
the peak in its kick profile; the relative dela y in merger times b e tween d ata sets differing in 
initial spins is consistent with the results of lCampanelli et~~al. (2006rf). All configurations 

show a marked "unkick" a fter the peak, with the exception of the equal-mass case, EQ H . 

[from 



Baker et al\ (|2007af ) - reproduced with permission of the AAS] 



calculate the merger waveforms accurately and robustly. Today there is consensus that 
the merger of equal-mass, nonspinning black holes produces a final Kerr black hole with 
spin a/rrif ~ 0.7, and that the amount of energy radiated in the form of gravitational 
waves, starting with the final few orbits and proceeding through the plunge, merger and 
ringdown, is £rad = M — rrif ~ 0.04M. Simulations are now being carried out for an 
increasing range of black-hole component masses and spins. Computational methods con- 
tinue to improve, resulting in longer simulations that start in the late inspiral regime and 
allow quantitative comparisons with analytic post-Newtonian methods. Applications of 
the resulting waveforms to gravitational wave data analysis have begun. The recoil or 
kick velocity that results when the black holes have unequal masses, or unequal or un- 
aligned spins, has been calculated for a variety of interesting cases, with key applications 
to astrophysics. Overall, the field of black-hole binary mergers is experiencing a true 
"golden age," with many new results coming out across a broad front. Stay tuned! 
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(Ames Research Center), and at the NASA Center for Computational Sciences (Goddard 
Space Flight Center). B.J.K. was supported by the NASA Postdoctoral Program at the 
Oak Ridge Associated Universities. S.T.M. was supported in part by the Leon A. Herreid 
Graduate Fellowship. 



12 



J. G. Baker et ai: Gravitational waves 



REFERENCES 

Abbott, B. et al. 2007 Search for gravitational- wave bursts in LIGO data from the fourth 

science run. arXiv:0704.0943 [gr-qc]. 
Alcubierre, M., Benger, W., Brugmann, B., Lanfermann, G., Nerger, L., Seidel, E. 

& Takahashi, R. 2001 3D Grazing Collision of Two Black Holes. Phys. Rev. Lett. 87, 

271103. 

Arnowitt, R., Deser, S. & Misner, C. W. 1962 The dynamics of general relativity. In 
Gravitation: An Introduction to Current Research (ed. L. Witten), pp. 227-265. New York: 
John Wiley. 

Baker, J., Campanelli, M., Pretorius, F. & Zlochower, Y. 2006a Comparisons of binary 
black hole merger waveforms. Class. Quantum Grav. 24, S25-S31. 

Baker, J. G., Boggs, W. D., Centrella, J., Kelly, B. J., McWilliams, S. T., Miller, 
M. C. & van Meter, J. R. 2007a Modeling kicks from the merger of non-precessing 
black-hole binaries. Astrophys. J. In press. 

Baker, J. G., Centrella, J., Choi, D.-L, Koppitz, M. & van Meter, J. 20066 Binary 
black hole merger dynamics and waveforms. Phys. Rev. D 73, 104002. 

Baker, J. G., Centrella, J., Choi, D.-L, Koppitz, M. & van Meter, J. 2006c Gravita- 
tional wave extraction from an inspiraling configuration of merging black holes. Phys. Rev. 
Lett. 96, 111102. 

Baker, J. G., Centrella, J., Choi, D.-L, Koppitz, M., van Meter, J. & Miller, M. C. 

2006a 7 Getting a kick out of numerical relativity. Astrophys. J. Lett. 653, L93. 
Baker, J. G., McWilliams, S. T., van Meter, J. R., Centrella, J., Choi, D.-L, Kelly, 

B. J. & Koppitz, M. 20076 Binary black hole late inspiral: Simulations for gravitational 

wave observations. Phys. Rev. D 75, 124024. 
Baker, J. G., van Meter, J. R., McWilliams, S. T., Centrella, J. fc Kelly, B. J . 2006e 

Consistency of post-Newtonian waveforms with numerical relativity, gr-qc/0612024 
Bardeen, J. M. & Petterson, J. A. 1975 The Lense-Thirring effect and accretion disks 

around Kerr black holes. Astrophys. J. 195, L65-L67. 
Belczynski, K., Kalogera, V., Rasio, F. A., Taam, R. E. & Bulik, T. 2007 On the rarity 

of double black hole binaries: Consequences for gravitational wave detection. Astrophys. J. 

662, 504-511. 

Blanchet, L. 2006 Gravitational radiation from post-Newtonian sources and inspiralling com- 
pact binaries. Living Rev. Rel. 9, 4. 

Bogdanovic, T., Reynolds, C. S. & Miller, M. C. 2007 Alignment of the spins of super- 
massive black holes prior to coalescence. Astrophys. J. 661, L147-L150. 

Brandt, S., Correll, R., Gomez, R., Huq, M. F., Laguna, P., Lehner, L., Marronetti, 
P., Matzner, R. A., Neilsen, D., Pullin, J., Schnetter, E., Shoemaker, D. & 
WlNlCOUR, J. 2000 Grazing collisions of black holes via the excision of singularities. Phys. 
Rev. Lett. 85, 5496-5499. 

Brugmann, B. 1999 Binary black hole mergers in 3D numerical relativity. Int. J. Mod. Phys. 
D 8, 85. 

Brugmann, B., Tichy, W. & Jansen, N. 2004 Numerical simulation of orbiting black holes. 

Phys. Rev. Lett. 92, 211101. 
Buonanno, A., Cook, G. B. & Pretorius, F. 2007a Inspiral, merger and ring-down of 

equal-mass black-hole binaries. Phys. Rev. D 75, 124018. 
Buonanno, A., Pan, Y., Baker, J. G., Centrella, J., Kelly, B. J., McWilliams, S. T. 

& VAN Meter, J. R. 20076 Toward faithful templates for non-spinning binary black holes 

using the effective-one-body approach. arXiv:0706.3732 [gr-qc]. 
Campanelli, M., Lousto, C. O., Marronetti, P. & Zlochower, Y. 2006a Accurate evo- 
lutions of orbiting black-hole binaries without excision. Phys. Rev. Lett. 96, 111101. 
Campanelli, M., Lousto, C. O. & Zlochower, Y. 20066 Last orbit of binary black holes. 

Phys. Rev. D 73, 061501(R). 
Campanelli, M., Lousto, C. O. & Zlochower, Y. 2006c Spin-orbit interactions in black-hole 

binaries. Phys. Rev. D 74, 084023. 
Campanelli, M., Lousto, C. O. & Zlochower, Y. 2006d Spinning-black-hole binaries: The 

orbital hang up. Phys. Rev. D 74, 041501(R). 



J. G. Baker et ai: Gravitational waves 



13 



Campanelli, M., Lousto, C. O., Zlochower, Y. & Merritt, D. 2007a Large merger recoils 
and spin flips from generic black-hole binaries. Astrophys. J. 659, L5-L8. 

Campanelli, M., Lousto, C. O., Zlochower, Y. k, Merritt, D. 20076 Maximum gravita- 
tional recoil. Phys. Rev. Lett. 98, 231102. 

Centrella, J. M. 2006 The final merger of comparable mass binary black holes. AIP Conf. 
Proc. 873, 70-79. 

Echeverria, F. 1989 Gravitational- wave measurements of mass and angular momentum of a 

black hole. Phys. Rev. D 40, 3194. 
Fitchett, M. J. 1983 The influence of gravitational wave momentum losses on the centre of 

mass motion of a Newtonian binary system. Mon. Not. R. Astron. Soc. 203, 1049-1062. 
Flanagan, E. E. & Hughes, S. A. 1998 Measuring gravitational waves from binary black hole 

coalescence: I. signal to noise for inspiral, merger, and ringdown. Phys. Rev. D 57, 4535. 
Fregeau, J. M., Larson, S. L., Miller, M. C, O'Shaughnessy, R. & Rasio, F. A. 2006 

Observing IMBH-IMBH binary coalescences via gravitational radiation. Astrophys. J. 646, 

L135. 

Gonzalez, J. A., Hannam, M. D., Sperhake, U., Brugmann, B. & Husa, S. 2007a Super- 
massive recoil velocities for binary black-hole mergers with antialigned spins. Phys. Rev. 
Lett. 98, 231101. 

Gonzalez, J. A., Sperhake, U., Brugmann, B., Hannam, M. & Husa, S. 20076 Maximum 
kick from nonspinning black-hole binary inspiral. Phys. Rev. Lett. 98, 091101. 

Hahn, S. G. & Lindquist, R. W. 1964 The two body problem in geometrodynamics. Ann. 
Phys. 29, 304-331. 

Hannam, M., Husa, S., Sperhake, U., Brugmann, B. & Gonzalez, J. A. 2007 Where 
post-Newtonian and numerical-relativity waveforms meet. arXiv:0706.1305 [gr-qc]. 

Herrmann, F., Hinder, I., Shoemaker, D., Laguna, P. fc Matzner, R. A. 2007 Gravita- 
tional recoil from spinning binary black hole mergers, gr-qc/0701143 

Herrmann, F., Shoemaker, D. & Laguna, P. 2006 Unequal-mass binary black hole inspirals 

Husa, S., Gonzalez, J. A., Hannam, M., Brugmann, B. & Sperhake, U. 2007a Reducing 
phase error in long numerical binary black hole evolutions with sixth order finite differenc- 
ing. larXiv:0706.0740l [gr-qc] . 

Husa, S., Hannam, M., Gonzalez, J. A., Sperhake, U. & Brugmann, B. 20076 Reducing 
eccentricity in black-hole binary evolutions with initial parameters from post-Newtonian 
inspiral. larXr7: 0706.0904 [gr-qc]. 

Kidder, L. E. 1995 Coalescing binary systems of compact objects to (post) 5//2 -Newtonian order. 

V. Spin effects. Phys. Rev. D 52, 821-847. 
Koppitz, M., Pollney, D., Reisswig, C, Rezzolla, L., Thornburg, J., Diener, P. & 

Schnetter, E. 2007 Recoil velocities from equal-mass binary-black-hole mergers. Phys. 

Rev. Lett. 99, 041102. 

Lang, R. N. & Hughes, S. A. 2006 Measuring coalescing massive binary black holes with 

gravitational waves: The impact of spin-induced precession. Phys. Rev. D 74, 122001. 
Leaver, E. W. 1986 An analytic representation for the quasi-normal modes of Kerr black holes. 

Proc. R. Soc. London Ser. A 402 (1823). 
van Meter, J. R., Baker, J. G., Koppitz, M. & Choi, D.-I. 2006 How to move a black hole 

without excision: gauge conditions for the numerical evolution of a moving puncture. Phys. 

Rev. D 73, 124011. 

Milosavljevic, M. & Merritt, D. 2003 Long term evolution of massive black hole binaries. 
Astrophys. J. 596, 860. 

Misner, C. W., Thorne, K. S. & Wheeler, J. A. 1973 Gravitation. San Francisco: W. H. 
Freeman. 

Pan, Y., Buonanno, A., Baker, J. G., Centrella, J., Kelly, B. J., McWilliams, S. T., 
Pretorius, F. & van Meter, J. R. 2007 A data-analysis driven comparison of analytic 
and numerical coalescing binary waveforms: non-spinning case. Submitted to Phys. Rev. 
D. 

Pfeiffer, H. P. et al. 2007 Reducing orbital eccentricity in binary black hole simulations. 

Class. Quantum Grav. 24, S59-S82. 
Pretorius, F. 2005 Evolution of binary black hole spacetimes. Phys. Rev. Lett. 95, 121101. 



14 



J. G. Baker et ai: Gravitational waves 



Schnittman, J. D. & BuONANNO, A. 2007 The distribution of recoil velocities from merging 
black holes. |astro-ph/0702641| 

Sesana, A., Haardt, F., Madau, P. & Volonteri, M. 2005 The contribution of inspiraling 
massive black hole binaries to the LISA data stream. Astrophys. J. 623, 23-30. 

Shaddock, D. A., Tinto, M., Estabrook, F. B. & Armstrong, J. W. 2003 Data combi- 
nations accounting for LISA spacecraft motion. Phys. Rev. D 68, 061303. 

Shoemaker, D. 2006 Private communication. 

Smarr, L. 1977 Spacetimes generated by computers: Black holes with gravitational radiation. 
Ann. N. Y. Acad. Sci. 302, 569-604. 

Smarr, L. 1979 Gauge conditions, radiation formulae and the two black hole collision. In Sources 
of Gravitational Radiation (ed. L. Smarr), p. 245. Cambridge, England: Cambridge Uni- 
versity Press. 

Smarr, L., Cadez, A., DeWitt, B. & Eppley, K. 1976 Collision of two black holes: Theo- 
retical framework. Phys. Rev. D 14 (10), 2443-2452. 

Tichy, W. & Marronet ti, P. 2007 Binary black hole mergers: Large kicks for generic spin 
orientations, gr-qc/0703075 

Vaishnav, B., Hinder, I., Herrmann, F. & Shoemaker, D. 2007 Matched filtering of nu- 
merical relativity templates of spinning binary black holes. arXiv:0705.3829 [gr-qc]. 



